Space filling in the World Trade Web: measures and null models 
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It is here proposed global and local, bynary and weighted measures of space filling, to detect the 
degree of stretching of a network within its embedding space. Measures are normalized indexes 
that vary from zero - a maximally shrunk network, to one -a maximally expanded network. We 
evaluated the indexes according to three different null models an propose an unbiased measure 
of network filling, suitable for confrontation of different (cross sectional) embedded networks or 
historical trends of the same network. We further develop a local measure of filling and thereby of 
vertex assortativity. Finally, we showed that in the WTW metric distances seem to have determined 
the links' allocation more than the present null models tend to have predicted and that this burden 
upon links has been changing in time. However, surprisingly, our analysis seems to indicate that 
during globalization the WTW expanded merely because of topology and contrary to expectations 
it actually shrunk. 



INTRODUCTION 



In recent years transport networks have prooved to be 
an insightful approach to some prominent and apparently 
unrelated scientific questions concerning the size and the 
form of recurrent structures in Nature ^0]. The ubiq- 
uitous allometric scaling law or the shape of some living 
and non-living systems have been demonstrated to be the 
result of efficient and optimized transport networks [5]- 
rTU] . More broadly, in the field of ecology, networks have 
been already approached as transport systems [TTJ . 
Indeed, the study of complex networks too has a long 
and successful record, that fostered progress in several 
different disciplines [T5lU7| . More recently, network the- 
ory has been applied to the study of transport networks 
in ecology [TBI EH] an d to the study of some cargo ship- 
ments and flights [20, 21 . However, so far, statistical me- 
chanics of networks has never been extensively applied to 
the study of transport networks. A transport network is 
a network embedded in an metric space, where vertexes 
can be positioned at fixed coordinates or free to move. 
In this paper we propose a method to study, by means of 
a statistical mechanics approach, to what extent an em- 
bedded network, with fixed vertexes, fills its embedding 
space. 



DEFINITIONS 

In a binary, directed and embedded network a given 
number of links L can be placed in between close or far 
apart vertexes. In the former case the network tends to 
be contracted and in the the latter tends to be stretched 
across its emebbedding space, or to fill it. Metric dis- 
tances in an embedded network are defined for every dyad 
(coulpe of nodes) and they can be orderly arranged from 
the smallest to the largest, for all the possible n{n — 1) 
combinations, and vice versa: = (d\, d2..-d n ( n -i)) 
where di = d 2 < d 3 = d 4 < ..d n ( n _i) and d} = 
(di,d 2 — rf„(„_i)) where di = d 2 > d 3 = d 4 > ..d n („_i). 
Links' placement can therefore span between two ex- 
tremes, the vector of the minimal distances or the vector 
of the maximal distances: min = d\ and max = 
dj . Summation over all the possible occupied dis- 
tances goes from 1 to L as this is the number of links 
defining the graph G* and, as it will be clearer later, it 
is also the expected number of links of the conjugated 
null model(s); thus, the constrain generating the graph 
ensemble G. It is thereby proposed global index of the 
space filling of an embedded network: 
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where a^- are the entries of the adjacency matrix. For 
any given number of links L this index can vary between 
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and 1, where the two extremes stand for a maximally 
shrink and a maximally expanded network. It is note- 
worthy that this index does not depend on the euclidean 
dimension of the embedding space, nor it depends on the 
topology of the network or the shape of the two (three) 
dimensional embedding object, whose boundaries are de- 
limited by the peripheral (outer) vertexes. It is further 
possible to derive a definition of network filling for the 
weighted matrix. If we intend weights as multiple links of 
unitary weight, whose statistical properties are described 
by a Bose-Einstein distribution [22], the two extremal 
configurations become: the whole weight W placed upon 
the closest dyad or the farthest. The weighted filling in- 
dex would thus be: 
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where: 
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w = y. 



Null Models 

The definition of network filling as expressed in equa- 
tion [T] assumes that links can be freely allocated in the 
network, that is, that there is a constant probability 
Pij = p of a linked to be formed between any dyad \l'dV 
]. Therefore, the conjugated null model is described by 
the following Hamiltonian: 



i?i(G) = 9L(G) =6^1 
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However, real networks are much more heterogeneous 
and featured by power law distributions of the main topo- 
logical properties. Scale free networks, for example, have 
a broad distribution of in and out degree |16j . Therefore, 
we want to measure the [T] in such a way that the real de- 
gree distribution in the network is preserved by the links 
rewiring. That is to say, we want that each vertex' con- 
tribution to the maximal or minimal link 's length is to be 
consistent with the observed topology. In this latter case 
the underlying null model is the directed configuration 
model, defined by the enforcement of both the in-degree 
and the out-degree sequences {fc-™} and {k° ut } as con- 
straints I23H251: 



H 2 {G) = WkTiG) + 9° ut k° ut {G)} 

i 

A null model tells us what are the trivial (maximally en- 
tropic) properties of the graph ensemble obtained by im- 
posing determined constrains into the hamiltonian. We 
can thus compare our observed properties to the excepted 



ones in order to set a level of triviality in the network 
for the detected property. That is to say: we can dis- 
entangled the property under investigation from other 
topological effects. In what follows every network, re- 
gardless of its topological features can be compared with 
one other. Moreover, a network can be detected along 
its evolution in time. This latter method was proposed 
to investigate the reciprocity in binary networks, where 
link density decisively affects the probability of mutual 
links [215] to form. The observed value can be confronted 
with the expected value according to a specified rewiring 
algorithm according to the following relation: 
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Hitherto, a varies from —1 to 1 and it scores whenever 
the observed filling coincides with the expected filling, 
according to the specified model. It is interesting that ac- 
cording the Hamiltonian defined in equetion[3]and to the 
one defined in equation [3j the conditional probability of a 
link to reciprocate another link is constant [23]. However, 
in a network where distances have a positive or negative 
effect upon vertex' connecting decisions, the conditional 
probability of links would arguably be affected. It is very 
likely that distances would possibly increment the sym- 
metry of the entries of the adjacency matrix. It is there- 
fore of our interest to require that the ensemble of our 
null model is to be symmetrical under link reversal per- 
mutation. One way to achieve this latter symmetry is to 
impose the reciprocal degree sequence to the ensemble: 

h 3 (g) = Y J [ e TkT{G) + er t K ut {G) + erkr{G)\ 
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Finally, it might be interesting to evaluate a third kind 
of null model, still symmetrical to transpose equivalence, 
but less restrictive compared to the former one (equation 
H). This model relaxes the constrains on every vertex' 
reciprocal degree, but preserves the reciprocity structure 
of the network across the ensemble globally: 

H 4 {G) = [ d TK l {G) + 9° ut K ut {G) + 9 r(G)] 
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This latter model is interesting because there are some 
networks, like the World Trade Web, that present a uni- 
form conditional probability, as it has been shown that 
the reciprocated degree is proportional to total degree for 
every vertex [25] and that this relations has bee steady 
over time. 
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Local filling and assortativity 

From the analysis of the weighted matrix it is further 
possible to derive a local network filling, inward and out- 
ward, and thereby a in and out assortativity. Assorta- 
tivity is a measure of the tendency of vertexes with de- 
termined common topological properties to connect with 
one other [27j. If we define a min /oc = £f /n 4(Ej w *i + 

J2 3 Wji) and max' oc = J2i d i(J2j w ij + Ej w ji)> where 
the summation over the distances neds at L/n as this lat- 
ter is the expected number of links according to the con- 
jugated null model, the Weigthed Random Graph [25] . 

We thereby proposed the following measure for the lo- 
cal filling: 



l out {i) = 
and a weighted version: 
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In figures [T] and [4] are plotted the observed assortativity 
and the in/out degree in the World Trade Web compared 
to the predicted assortativity according to the Weighted 
Configuration Model [25] . 



RESULTS 

We have evaluated the network filling for the (binary 
and weighted) World Trade Web (data source: UN COM- 
TRADE). The WTW, exhibits a constant growing trend 
in volume and connectance, from 86 vertexes and a 0.39 
link density in 1950 to 187 and 0.57 in 2000. Interest- 
ingly, the average distance and variance have been pretty 
stable all over the time, going from 7516 and 2.3* 10 7 Km 
in 1950, to 7550 and 1.9* 10 7 Km in 2000. However, net- 
work filling was not steady and shows a growing trend in 
the time that goes from the late 1950s to the late 1970s 
and a following decreasing trend (figure [5]) . Surprisingly 
it seems that during globalization the network actually 
shrunk. Nevertheless, it should be noted that the spatial 
expansion of the WTW seems to coincide with the sudden 
and unexplained leap in reciprocity occurred in the WTW 
during the lag of time that goes from the late 1980s on- 
ward [25]. This is noteworthy because, in spite of the fact 
that the directed configuration model well explains most 
of the topological properties of the WTW, reciprocity re- 
mains unexplained [30] , Furthermore, as was previously 
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FIG. 1: Out — assortativity and out — degree, in the 
binary WTW (Year 2000), observed (Blue)and predicted 
with DCM(Red). 
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FIG. 2: Local out — filling and out — degree, in the binary 
WTW (Year 2000), observed (Blue) and predicted with 
DCM(Red). 



highlighted, reciprocity arguably affects, or it is affected 
by metric distances of dyads. Therefore, it might be pos- 
sible that including the reciprocity structure locally (H3) 
or globally (H4) is to be enlightening in some way. In fig- 
ure [5] the four null models collapse into two, indicating 
on the one hand that null models based on one global 
parameter {Hi and intend to behave very similarly; on 
the other hand that the local reciprocity structure (-H4) 
does not add significant information to the ensemble in 
order to explain the space filling compared to the degree 
sequence {Hi). However, surprisingly, both these two lat- 
ter null models indicate that the space filling, discounted 
of the topological effects, actually increased during the 
previously mentioned lag of time, as depicted in figure 
g] That is to say, after late 1980s, the WTW expanded 
and this phenomenon was arguably caused by a change 
in the topology rather than other factors affecting trades. 
Interestingly, the weighted filling of the WTW seems to 
display the same opposite trend: a sudden collapse in 
the weighted filling is detectable during the period late 
1970s early 1990s (figure [7]). Hence, while the binary 
structure of the network expanded, but merely because 
of its topology, the weighted shrunk. This apparent di- 
chotomy is interesting in the light of the divergent results 
of the binary an weighted (out) assortativity. Despite lo- 
cally the weighted and binary version of the graph behave 
similarly (figures [2] and [4), the weighted graph is assor- 
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FIG. 3: Out — assortativity and out — strenght, in the 
weighted WTW (Year 2000), observed (Blue) and pre- 
dicted with WCM (Red). 
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FIG. 6: Global <f> evaluated according to the Reciprocated 
DCM {Hi) 




FIG. 4: Local out — filling and out — strenght, in the 
weighted WTW (Year 2000), observed (Blue) and pre- 
dicted with WCM (Red). 



FIG. 7: Weighted filling in the WTW: observed 
(blue), WCM (yellow), WRG (red) 



tative whereas the binary representation is disassortative 
(figures [l] and [3| . That is, though far apart (peripheral 
in the metric space) vertexes tend to be connected with 
hubs (central vertexes in the topological space), vertexes 
hilly connected tend to place their weights prominently 
with closed (neighbours) vertexes. 
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FIG. 5: Expected and observed space filling in the WTW: 
Hi (red), H2 (orange), H3 (yellow), Hi (green) and ob- 
served /* (blue) 

A further interesting result of the analysis is that all 
the four null models evaluated tend to over estimate the 
value of network filling for the WTW. The WTW should 
be a shrunk network, where distances pose a significant 
burden over trade decisions. There is a long standing 



debate in literature over the role of distances in trades 
[5TI [33] and more recently, Fagiolo tested the effect of 
distances on the topology of the WTW [33] . Our results 
seem to confirm previous findings, in this respect. Nev- 
ertheless, it must be noted that locally both the WCM 
and the DCM, though still overestimating the space fill- 
ing, are more accurate, signalling once again that it is 
the global structure of the WTW that has a paramount 
effect on the degree of space filling of the network. We 
have also investigated the effect of distances on different 
scales of the network. In Table |T] it is shown that the 
European transportation network (data source: EURO- 
STAT) is more spatially expanded than the world wide 
network, whilst on a national scale, for the case of Italy 
(data source: ISTAT), the network seems to be signif- 
icantly shrunk {<f> = —0.142). This picture gives the 
impression of a (binary) network where distances count 
more on a local scale than a global scale, but this result 
might be significantly affected by the sources of data and 
more research is needed. 

TABLE I: Space filling in transport networks 





WTW (2000) 


EU 15 (2000) 


Italy (2010) 


Link density 


0.578 


0.819 


0.844 


f 


0.578 


0.458 


0.377 


<t> (CM) 


0.380 


0.134 


-0.142 


r 


0.228 


0.379 


0.328 
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CONCLUSIONS 

We have here proposed a possible measures to detect 
the space filling, that is to say the extent of stretching, of 
a space embedded network and confronted this with four 
different null models, in order to develop an improved 
version of the index. Despite the rough version of the in- 
dex seems to indicate that during globalization the net- 
work, contrary to expectations, contracted, the improved 
version, discounting topological effects, suggest that the 
WTW expanded significantly after the late 1980s. Our 
analysis further shows that reciprocity is undeniably con- 
nected to the space filling, though it does not explain 
it. It is arguably possible that the same unknown phe- 
nomenon that affected globally the reciprocal structure 
of the WTW also caused its spatial expansion. This hy- 
pothesis is further confirmed by the local analysis of the 
space filling, that is fairly well described and explained, 
contrary to the global analysis, by the considered null 
models. We have to deduce, therefore, that this unknown 
phenomenon was global and structural, yet more research 
is needed to detect it. 
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